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Null test of the cosmic curvature using H{z) and supernovae data 
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We introduce a model-independent approach to the null test of the cosmic curvature which is 
geometrically related to the Hubble parameter H{z) and luminosity distance dciz). Combining the 
independent observations of H{z) and driz), we use the model-independent smoothing technique, 
Gaussian processes, to reconstruct them and determine the cosmic curvature 0^^ in the null test 
relation. The null test is totally geometrical and does not assume any cosmological model. We show 
that the cosmic curvature = 0 is consistent with current observational data sets, falling within 
the la limit. To demonstrate the effect on the precision of the null test, we produce a series of 
simulated data of the models with different Future observations in better quality can provide 

a greater improvement to constrain or refute the fiat universe with = 0. 

PACS numbers: 98.80.-k, 98.80.Es, 98.80Jk 


I. INTRODUCTION 


Whether the space of our Universe is open, flat, or 
closed is one of the most fundamental problems in mod¬ 
ern cosmology. The spatial topology of the Universe is 
also closely related to other important problems such as 
the evolution of the Universe, the property of the dark 
energy, etc. The effect of allowing nonzero curvature on 
dark energy models has been studied in many papers, see, 
for example, Refs. And the detection of a signifi¬ 

cant deviation from = 0 would have profound con¬ 
sequences for inflation models and fundamental phwics. 
A lot of attention has been attracted to this issue [§-Q. 
The AGDM model is consistent with all of the data and 
a flat universe is preferred even in recent Planck 2015 
results 0. 

However, almost all of the papers studying and con¬ 
straining the cosmic curvature assume some specific mod¬ 
els for dark energy such as the equation of state w{z). 
These are all model-dependent and indirect methods. 
Note the fact that there is some degeneracy between 
the spatial curvature and the equation of state of dark 
energy in these studies. It would be better to detect 
the space curvature of the Universe by geometrical and 
model-independent methods. In this paper, we use a 
geometrical relation among the cosmic curvature 
Hubble parameter H{z), and luminosity distance dL{z), 
by combining the Hubble rate and luminosity distance, 
we are able to directly determine and test whether the 
cosmic curvature deviates from 0 [E 113) Ell- For this 
goal, we should focus on two independent observations 
that directly give H[z) and dL{z), respectively. For H{z) 
data, it can be derived from differential ages of galaxies 
[“cosmic chronometer (GG) ”] and from the radial baryon 


acoustic oscillation (BAG) scale in the galaxy distribu¬ 
tion. As for dL{z), we use the SNIa Union 2.1 data sets. 
We use the model-independent method Gaussian pro¬ 
cesses (GP) for smoothing the observational data. The 
advantage of these methods is that they are all model- 
independent; hence, we need not assume any models in¬ 
volving dark energy and gravity theory. They are purely 
geometrical and are constrained directly by observational 
data. However, the precision of the null test is also lim¬ 
ited by the quality of observational data. If we want to 
detect a tiny cosmic curvature more precisely, a better 
quality of the observational data sets are also required, 
which will be also discussed in this paper. 

This paper is organized as follows. In Sec. im we intro¬ 
duce the theoretical method for the null test of the cosmic 
curvature. In Sec. imi we first give a brief introduction 
of Gaussian processes, and then apply GP method to the 
null test of the cosmic curvature using two independent 
data sets: CC-I-BAO, Union 2.1. Further it is followed 
by a series of simulated data tests. We give discussions 
and conclusions in Sec. El 


II. THEORETICAL METHOD 


In a FLRW universe, the luminosity distance can 
be expressed as 


di = 


c(l -I- z) 

Ho^ 


sinh( 



dz ^ 

Wr 


( 1 ) 


where E{z) = H{z)/Hq, = —Kc^/{aQHQ)^, and 
K = -1-1,—1,0 corresponds to a closed, open, and flat 
universe, respectively. 

Differentiating Eq. o and writing D{z) = (ido/c)(l -I- 
z)~^dL{z) as the normalized comoving distance, we can 
obtain 


n 


(0) 

K 


E^{z)D'^{z) - 1 
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We can see that the cosmic curvature can be directly 
determined by using the Hubble parameter and luminos¬ 
ity distance from Eq. ([5|). Thus, the null test of is 
straightforward. Note the fact that D{0) = 0 will bring 
about a singularity at 2 ; = 0. Therefore for more succinct, 
we transform Eq. m to 


E{z)D'{z) + 1 


E{z)D'{z) - 1. 


( 3 ) 


Since the left-hand side (Ihs) of Eq. ([3]) is a nonzero when 
2 ; 0 if is nonvanishing, the null test of cosmic 


curvature is equivalent to the null test of the whole 
Ihs of Eq. ([3]) . If we define 


Ok{z) 


E{z)D'{z) -f 1’ 


then a flat universe implies 


Ok{z) =E{z)D'{z) -1 = 0 


( 4 ) 

( 5 ) 


is always true at any redshift. A deviation from it will 
indicate a signal of nonvanishing cosmic curvature. Note 
that the theoretical value of Ok{z) is always zero at 2 ; = 
0. So the null test is just to check whether there exists 
signal of nonvanishing Ok{z) at nonzero redshifts. 

Thus we should use current observational data sets to 
reconstruct E{z) and D'{z), independently, and combine 
these two reconstructions to test whether the relation 
E{z)D'{z) — 1=0 holds at arbitrary redshifts. We want 
to stress here that the null test is totally geometrical and 
cosmological model independent, the reconstructions of 
E{z) and D'{z) are directly derived from observational 
data sets. 


III. NULL TEST USING H{z) AND 
SUPERNOVAE DATA 

Given some observational data sets, it is crucial to 
use a mo del-independent method to reconstruct E{z), 
D{z), and its derivative D'{z). There are many dif¬ 
ferent methodologies to reconstruct functions from the 
data. For a brief analysis, see [ll]. Since we should 
use a nonparametric approach to smooth the data and 
to reconstruct the derivative, the so-called Gaussian pro¬ 
cesses [T3l - [l^ are very suitable for our purpose. 


A. Gaussian processes 


observational data. At each point 2 :, the reconstructed 
function f{z) is also a Gaussian distribution with a mean 
value and Gaussian error. The functions at different 
points z and z are related by a covariance function A:(z, z), 
which only depends on a set of hyperparameters i and 
CT/. Here, i gives a measure of the coherence length of 
the correlation in the x direction and cr/ denotes the over¬ 
all amplitude of the correlation in the y direction. Both 
of them will be optimized by GP with the observational 
data sets. In contrast to actual parameters, GP does not 
specify the form of the reconstructed function. Instead, 
it characterizes the typical changes of the function. The 
detailed anaHsis and description of the GP method can 
be found in [l^, [l^ . 


B. Hubble rate data and Union 2.1 

Following [l3, [III, we proceed to an analysis based 
on observational Hubble data compiled from several 
sources, independent of SNela. We combine measure¬ 
ments of H{z) obtained with two methods. One is cos¬ 
mic chronometers, which are mainly passively evolving 
galaxies. There are 21 data points compiled by Moresco 
et al. [H, [2II . The other is radial baryon acoustic oscil¬ 
lations from galaxy clustering in redshift surveys, which 
gives seven data points of Hubble parameters from dif¬ 
ferent experiments [ItI - I^ . We summarize the total 28 
data points in Table HI 

We normalize H{z) using E[q = 70km/(s Mpc); thus, 
we get the observational data points of E{z). Then we 
can use GP method to reconstruct E{z). Note that iJo is 
just a normalization factor, whose value will not influence 
our null test Eq. ([5]). 

To reconstruct U(z), we use SNela Union 2.1 data 
sets [HI, which contain 580 SNela data. We transform 
the distance modulus m — M given in the data set to D 
using 


m — M + 5 log 



-25 = 5 log [{l + z)D]. 


( 6 ) 


For consistency, here we also use Hq = 70km/(s Mpc) to 
normalize dL{z). Obtaining these 580 observational data 
points of D{z), we can also use GP method to reconstruct 
D{z) and its derivative D'{z). Finally, we combine the 
reconstructions of E{z) and D'{z) and apply them to the 
null test of Ok{z) in Eq. ([H). We stress again that the 
null test is model independent, so we need not assume 
any cosmological model, and the two data sets of H{z) 
and supernovae are also independent of each other. 


The Gaussian processes allow one to reconstruct a 
function from data without assuming a parametrization 
for it. Here, we use the Gaussian processes in Python 
(GaPP) [ 1 ^. This GP code has been used in various 
papers for different studies [Tbilisi. The distribution 
over functions provided by GP is suitable to describe the 


C. Null test 

Having obtained the total 28 data points of E{z) and 
580 points of D{z), we now use the GP method to recon¬ 
struct them, respectively. 
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FIG. 1: Gaussian precess reconstruction of E{z) (left) from GC+BAO, D{z) (middle), and D'{z) (right) from Union 2.1. The 
shaded blue regions are the 68% and 95% C.L. of the reconstruction. The flat ACDM model (red line) with flrno = 0.27 is also 
shown. 


Index 

z 

Hiz) 

Refs. 

1 

0.090 

69± 12 

[25] 

2 

0.170 

83 ±8 

[2^ 

3 

0.179 

75 ±4 

[25] 

4 

0.199 

75 ±5 

[25] 

5 

0.240 

79.69 ± 2.32 

[27] 

6 

0.270 

77± 14 

m 

7 

0.350 

82.1 ±4.9 

[2^ 

8 

0.352 

83± 14 

m 

9 

0.400 

95 ±17 

m 

10 

0.430 

86.45 ±3.27 

[27] 

11 

0.440 

82.6 ±7.8 

[29] 

12 

0.480 

97 ±62 

m 

13 

0.570 

92.4 ±4.5 

[30] 

14 

0.593 

104 ± 13 

[2^ 

15 

0.600 

87.9 ±6.1 

[29] 

16 

0.680 

92 ±8 

[2^ 

17 

0.730 

97.3 ± 7 

[29] 

18 

0.781 

105 ± 12 

[2^ 

19 

0.875 

125 ± 17 

[2^ 

20 

0.880 

90 ±40 

[2^ 

21 

0.900 

117 ±23 

[2^ 

22 

1.037 

154 ± 20 

[2^ 

23 

1.300 

168 ± 17 

[2^ 

24 

1.363 

160 ± 33.6 

[26] 

25 

1.430 

177 ±18 

[25] 

26 

1.530 

140 ± 14 

[25] 

27 

1.750 

202 ± 40 

[25] 

28 

1.965 

186.5 ± 50.4 

m 
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FIG. 2: Reconstruction of Ok{z) from CC+BAO and Union 
2.1. The shaded blue regions are the 68% and 95% C.L. of the 
reconstruction. The red line corresponds to the flat universe 


n 


( 0 ) 


= 0 . 


higher redshifts, the errors become large due to the poor 
quality data in that region. Using the reconstructions of 
E{z) and D'{z), we apply Monte Carlo sampling to de¬ 
termine the Ok(z) in Eq. ([S]) at each point z, which we 
want to reconstruct. In Fig. [2l it is shown that the recon¬ 
structed Ok(z) is consistent with the vanishing cosmic 
curvature, falling in the la limit. It tells us that there is 
no significant signal to indicate the deviation of the cos¬ 
mic curvature from 0 at the current observational 
data [H{z) and supernovae] level. In addition, let us 
mention that the mean value of the cosmic curvature is 
negative in the high redshift region, which is also consis¬ 
tent with the results from model-dependent constraints 
in the literature. 


TABLE I: H{z) measurements from different surveys using 
passively evolving galaxies and radial BAO. 


D. Mock data 


We can see from Fig. [Tjthat all of the reconstructions 
of E{z), D{z), and D'{z) are consistent very well with 
the flat ACDM model, which we assume for comparison. 
The dashed blue line is the mean of the reconstruction 
and the shaded blue regions are the 68% and 95% confi¬ 
dence level (C.L.) of the reconstruction. As expected, at 


To demonstrate how a large number of data with what 
accuracy of the error will affect our null test when re¬ 
constructing E[z) and D[z), we firstly simulate a data 
set of 128 points for E{z). Adopting the methodol¬ 
ogy in we draw the error from a Gaussian dis¬ 

tribution: ue ^ Af{a,e) with a = (cr+ -I- a-)l2 and 
e = {a+ — cr_)/4, where cr+ and cr_ are the two straight 
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FIG. 3: Gaussian process reconstruction of E{z), D{z) (top), D'{z), and reconstruction of OK{z){bottom) obtained from a 
mock data set of E{z) and future DES, assuming the concordance model with = 0 (red line). The dashed blue line is the 
mean of the reconstruction, and the shaded blue regions are the 68% and 95% C.L. of the reconstruction, respectively. 




lines that bound the uncertainties cr^z) of the observa¬ 
tional E(z) data from above and below, respectively. 
Then E{z)sim is sampled from the Gaussian distribution 
E{z)sim ^ Af{E{z)fid, (Je), where E{z)fid is the theoret¬ 
ical value from the fiducial model. 

As for simulated D{z) data, we create mock data sets 
of future SNela according to the Dark Energy Survey 
(DES) [1^. The DES is expected to obtain high qual¬ 
ity light curves for about 4000 SNelas from z = 0.05 
to z = 1.2. From Table 14 in [s^, we can calcu¬ 
late the errors of D, ud, and the corresponding num¬ 
bers of SNela for each redshift bin. At every redshift 
point z, D[z)sim is sampled from the normal distribu¬ 
tion D{z)sim Af{D{z)fid,crD)- 

Having obtained the simulated data sets of E{z) and 
D{z), we use GP method to reconstruct E{z), D{z), and 
D'{z). Then we combine the reconstructions of E{z) 
and D'(z), using Monte Carlo sampling to determine the 
Ok{z) and check whether GP can recover its theoretical 
value and distinguish it from = 0. 

We now simulate the data points for three different 
fiducial models: concordance model, namely, ACDM 
model with = 0 and = 0.3; two models with non¬ 
vanishing cosmic curvature, = ±0.16 and = 0.3. 
We want to check whether the GP method can detect or 
recover all of them and distinguish them from each other. 
The results are shown in Figs. [51[S1 respectively. 

We can see from Fig.[3]that E(z), D(z), and D'{z) are 
all reconstructed very well from the mock data sets as¬ 
suming the concordance model = 0. And the recon¬ 


structed Ok{z) is also consistent with the concordance 
model nicely. As expected, at higher redshifts, the errors 
become large due to the poor quality data in that region. 
Furthermore, we see from Figs. H] and [5] that the recon¬ 
structions of model = ±0.16 also recover the fiducial 
model very well, falling in the Icr limit and obviously de¬ 
viating from the concordance model = 0 at 95% 
G.L. The large errors and not so good reconstructions at 
high redshifts (z > 1.0) are due to the poor quality data 
in that region for there are no simulated data of D{z) at 
z > 1.2. Anyway, we can claim that with these quantities 
and errors of the observational data for E{z) and D(z), 
the GP method at least has the ability to detect the cos¬ 
mic curvature > 0.16 and to rule out = 0 at 2cr 
G.L. 

However, as we have pointed out that the theoretical 
value of Ok{z) is always zero at z = 0. So the theoreti¬ 
cal values of Ok{z) for these models with different 

deviates from the model with D® = 0 tiny at low red¬ 
shifts. As z increases, the difference becomes large. As a 
result, we can see from Figs. 0] and [5] that it is very hard 
to rule out D® = 0 at low redshifts. Maybe it can work 
out at a middle redshift (0.6 < z < 1.0) as Figs. [Hand[5] 
show, but it is not always helpful when D® is smaller. 
So, if we want to detect a model with D® < 0.16 or even 

smaller and can rule out D® = 0 at 95% C.L., a larger 
number and higher quality data sets are required. 

Figure ini shows the reconstructed Ok{z) also using the 
mock data sets of E{z) and D{z) but with a quarter of 
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FIG. 4: Gaussian process reconstruction of E{z), D{z) {top), D'{z), and reconstruction of OK{z){bottom) obtained from a 
mock data set of E{z) and future DES, assuming the fiducial model with = 0.16 (green line). The dashed blue line is the 
mean of the reconstruction, and the shaded blue regions are the 68% and 95% G.L. of the reconstruction, respectively. The 
concordance model = 0 is also shown (red line). 





FIG. 5: Gaussian process reconstruction of E{z), D{z) {top), D'{z), and reconstruction of OK{z){bottom) obtained from a 
mock data set of E{z) and future DES, assuming the fiducial model with = —0.16 (green line). The dashed blue line is 
the mean of the reconstruction, and the shaded blue regions are the 68% and 95% C.L. of the reconstruction, respectively. The 
concordance model = 0 is also shown (red line). 
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FIG. 6: Reconstructions of Ok{z) for models with different the errors are a quarter of the original errors of the mock 
data sets, (a) = 0.05, (b) = —0.05, (c) = 0.09, (d) = —0.09, (e) = 0.16, (f) = —0.16, The shaded 

blue regions are the 68% and 95% C.L. for the reconstruction. 


the errors. The reconstructions are more nice. And the 
null test is more precise; it can detect the model with 
n® = 0.05 or even smaller, ruling out = 0 at 95% 

C.L. 


IV. DISCUSSIONS AND CONCLUSIONS 

In this paper, we have introduced a nonparametric ap¬ 
proach to making a null test of the cosmic curvature. Us¬ 
ing the Gaussian process method, we reconstructed the 
Hubble rate and distance-redshift relation \E{z), D{z), 


and D'{z)], independently. In the reconstruction, we 
needed not to assume any cosmological model. By com¬ 
bining the reconstructions of E(z) and D'{z), we can 
determine Ok{z)^ which is related to the cosmic curva¬ 
ture We have shown that 0^^ = 0 is consistent 

with current data sets (CC-I-BAO, Union 2.1), falling 
within the Icr limit, although the mean of the recon¬ 
structed curvature Ok is negative in the high redshift 
region, which is also consistent with the results from the 
model-dependent constraints in the literature. In addi¬ 
tion, note that the mean of the reconstructed curvature 
Ok in Fig. [3] is positive in the high redshift region for 
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the fiducial flat ACDM model. In this sense, our result 
shown in Fig 12] indicates that there is a little possibility 
for a closed universe. 

To demonstrate how a large number of data sets with 
different accuracies and errors will affect our null test, 
we create mock data sets of E(z) and p(z) using the 
methodology proposed in Refs. 32| and [^. We found 
that with current quality of simulated data, the GP 
method can recover and distinguish models with a differ¬ 
ent cosmic curvature with an order of 10“^ from 17^^ = 0. 


However, if we want to detect even smaller we 

should decrease the uncertainties of the data sets. Mak¬ 
ing the errors a quarter of the mock data, it can rule out 
= 0 from the cosmic curvature with an order of 10“^ 
at 95% C.L. 

Based on a one-parameter extension to the six- 
parameter ACDM model, CMB data d, HI, [l^ provide 


a strong constraint on . This constraint can be im¬ 
proved dramatically by adding BAO data that break the 
geometric degeneracy between and Hq. However, 
the model assumption and priors on model parameters 
may bias estimates of the cosmic curvature. Here, we 
proposed a model-independent method for reconstruct¬ 


ing the function Ok{z) with redshift, which character¬ 
izes a deviation from a flat universe. As we can see from 
Eq. (I5|), Ok{z) is completely and directly determined 
by the Hubble parameter and luminosity distance. If 
we can reconstruct the Hubble parameter and luminos¬ 
ity distance, the reconstruction of the cosmic curvature is 
straightforward. GP is suitable for our purpose because 
it can smooth data and reconstruct a function model in¬ 
dependently. Therefore, we can test the cosmic curvature 
neither assuming models nor imposing priors. Although 
current low-redshift data put weaker constraints on the 
cosmic curvature than GMB data, this method provides 
a new cross-check of the cosmic curvature using future 
data from large-scale structure measurements. 
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